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Abstract. Recent progresses in the study of jet modification in hot medium and their consequences in 
high-energy heavy-ion collisions are reviewed. In particular, I will discuss energy loss for propagating 
heavy quarks and the resulting modified fragmentation function. Medium modification of the parton frag- 
mentation function due to quark recombination are formulated within finite temperature field theory and 
their implication on the search for deconfined quark-gluon plasma is also discussed. 
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1 Introduction 

Jet quenching in high-energy nuclear collisions has been 
proposed as a good probe of the hot and dense medium 
formed in ultra-relativistic heavy-ion collisions. The 
quenching of an energetic parton is caused by multiple 
scattering and induced parton energy loss during its prop- 
agation through the hot QCD medium. Recent theoreti- 
cal estimates |3||3]Elini|7| all show that the effective par- 
ton energy loss is proportional to the gluon density of the 
medium. Therefore measurements of the parton energy 
loss will enable one to extract the initial gluon density of 
the produced hot medium. Because of color confinement in 
the vacuum, one can never separate hadrons fragmenting 
from the leading parton and particles materializing from 
the radiated gluons. The total energy in the conventionally 
defined jet cone in principle should not change due to in- 
duced radiation, assuming that most of the energy carried 
by radiative gluons remains inside the jet cone [5]. Addi- 
tional rescattering of the emitted gluon with the medium 
could broaden the jet cone significantly, thus reducing the 
energy in a fixed cone. However, fiuctuation of the un- 
derlying background in high-energy heavy-ion collisions 
makes it very difficult, if not impossible, to determine the 
energy of a jet on an event-by-event base with sufficient 
precision to discern a finite energy loss of the order of 10 
GeV. Since high-p^ hadrons in hadronic and nuclear colli- 
sions come from fragmentation of high-p^ jets, energy loss 
naturally leads to suppression of high-p-r hadron spectra 



El. 

Since parton energy loss effectively slows down the 
leading parton in a jet, a direct manifestation of jet quench- 
ing is the modification of the jet fragmentation function, 
Da^ft,(z, /i^), which can be measured directly in events 
in which one can identify the jet via a companion particle 
like a direct photon or a triggered high px hadron. This 
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modification can be directly translated into the energy loss 
of the leading parton. Since inclusive hadron spectra are a 
convolution of the jet production cross section and the jet 
fragmentation function in pQCD, the suppression of in- 
clusive high-pT hadron spectra is a direct consequence of 
the medium modification of the jet fragmentation function 
caused by parton energy loss. 

Strong suppression of high transverse momentum hadron 
spectra is indeed observed by experiments TD^, 1 1 at the 
Relativistic Heavy-Ion Collider (RHIC) at the Brookhaven 
National Laboratory (BNL), indicating large parton en- 
ergy loss in a medium with large initial gluon density. 
Shown in Fig.^are the nuclear modification factors RAAiPx) 
for single hadron spectra as a function of the number of 
participant nucleons. The theoretical results are obtained 
from a LO pQCD parton model calculation incor- 
porating modified parton fragmentation functions due to 
parton energy loss. 
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The initial gluon density in the most central Au + Au 
collisions at ^/s = 200 GeV was fixed by fitting the data. 
The centrality dependence shown is the consequence of the 
above parton energy loss, assuming the initial gluon den- 
sity is proportional to the measured hadron multiplicity 
which in turn is approximately proportional to the num- 
ber of participant nucleons. Such a calculation also has a 
definite prediction of the energy dependence of the hadron 
spectra suppression |I3| that agrees well with the current 
collection of data at different energies [T^ . 

Since the parton energy loss [Eq. JQ)] depends on the 
path length of the jet propagation which in turn depends 
on the azimuthal angle with respect to the reaction plane 
in non-central collisions, the parton energy loss and modi- 
fied fragmentation functions naturally lead to an azimuthal 
angle dependence of the hadron spectra suppression jl7l 
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Fig. 1. The centrality dep endence of the measured single in- + (lo*^'' curves) mid p + p (upper curves) colhsions as 

elusive hadron suppression at high-pT as compared to compared to the STAR|2T] data, 

theoretical calculation with parton energy loss 
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Fig. 2. Azimuthal anisotropy in Au + Au collisions as com- 
pared to the STAR .2U. 4-particle cumulant result. 

I18j . Indeed, the azimuthal angle distributions of high px 
hadrons were found to display large anisotropy with re- 
spect to the reaction planes of non-central Au + Au 
collisions. As shown in Fig. |S1 the observed azimuthal 
anisotropy, characterized by the second coefficient of the 
Fourier transformation V2, can also be described well by 
the same pQCD parton model calculation. 

The most striking measurement that is a direct mani- 
fest of jet quenching is the observed disappearance of the 
back-side high-p^ two-hadron correlation in azimuthal an- 
gle J which is characteristic of high-pT- back-to-back 
jets m. p -\- p collisions. Shown in Fig.|Sl are back-side two 
high Pt hadron correlations in Au + Au collisions with 
different centralities as compared to the same correlation 
in p + p collisions. Though the back-side correlations in 
peripheral Au + Au collisions remain the same as in p + p 



collisions, the peak gradually decreases and finally disap- 
pears in the most central Au + Au collisions. The curves 
are again from the same LO pQCD parton model calcu- 
lation with parton energy loss. It describes both the mag- 
nitude of the suppression and the centrality dependence 
quite well. 

Combining the above measurements of three different 
effects of parton energy loss and compare with the mea- 
sured jet quenching in deeply inelastic e + A collisions, 
one can conclude that the initial gluon density reached in 
central Au + Au collisions at y/s = 200 GeV is about 30 
times higher than in a cold Am nuclei [Hll^ . assuming the 
theoretical result that parton energy loss is proportional 
to the gluon density of the medium. 

In this talk, I will review recent progresses we have 
made in the study of modified jet fragmentation functions 
in medium and their applications to high px hadron spec- 
tra and correlations in heavy-ion collisions. In particular, 
I will review heavy quark energy loss and modified frag- 
mentation functions. Its unique features due to the heavy 
quark mass can help us to characterize the partonic nature 
of the observed jet quenching. I will also discuss dihadron 
fragmentation functions since it addresses hadron corre- 
lation within a jet and how it will also be modified by 
the parton energy loss. Finally, I will discuss the formu- 
lation of modified fragmentation functions due to quark 
recombination during the the hadronization of the jet in 
a thermal medium. The implication on the search for de- 
confined quark-gluon plasma will also be discussed. 
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2 Modified Heavy Quark Fragmentation 
Function 

In the study of parton energy loss, formation time for the 
radiated gluon plays an essential role. Because of the LPM 
interference, gluons with formation time longer than the 
length of a finite medium or the mean-free path in an in- 
finitely large medium will be suppressed. Such formation 
time is only relative to the propagation of the leading par- 
ton. Therefore, gluon radiation from a heavy quark is nor- 
mally shorter than that from a light quark because of the 
smaller velocity of the heavy quark. For radiated gluons 
with transverse momentum £± and z fractional momen- 
tum, the formation time is j23| 



2z{l-z)E 
~ 4 + (l-z)2M2^ 



(2) 



where M is the quark mass. Therefore, one should expect 
the LPM effect to be significantly reduced for intermediate 
energy heavy quarks. In addition, the heavy quark mass 
also suppresses gluon radiation amplitude at small angles 
121] relative to that off a light quark. Both mass effects 
will lead to a reduced heavy quark energy loss compared 
to that of a light quark. The most significant consequence 
of the reduced formation time due to heavy quark mass is 
the change of length dependence of the quark energy loss. 
The non-Abelian LPM effect due to suppression of gluon 
radiation with long formation time leads to a quadratic 
length dependence of the total energy loss. For a slow 
heavy quark, however, the dependence will become linear 
because of the short formation time and absence of the 
LPM interference. Shown in Fig. 0] are the numerical re- 
sults of the nuclear size Ra dependence of charm quark 
fractional energy loss in DIS off a cold nucleus, rescaled by 
C{Q'^)CAa'^{Q^)/Nc, for different values of xb and Q'^. 
One can clearly see that the Ra dependence is quadratic 
for large values of or small xb (large initial quark en- 
ergy) when the mass of the quark is negligible. The depen- 
dence becomes almost linear for small or large x b ■ The 
charm quark mass is set at M — 1.5 GeV in the numerical 
calculation. 

One can similarly calculate the nuclear modification 
of the heavy quark fragmentation as in the case of a light 
quark j22j . Shown in Fig. |S1 is the ratio of the modified 
charm quark fragmentation into D mesons to the vacuum 
fragmentation function. One can see that the modification 
due to the parton rescattering and induced gluon radia- 
tion in a nucleus for heavy quarks is quite different from 
light quarks |22ll25j . This is mainly caused by the form of 
the heavy quark fragmentation function in vacuum which 
peaks at large z. Because of the multiple scattering and 
induced gluon radiation, the position of the peak of the 
modified fragmentation function is effectively shifted to a 
smaller value of z. As a consequence, the heavy quark frag- 
mentation function remains unchanged, or even slightly 
enhanced for a large range of fractional momentum z. The 
modification only becomes significant and the fragmenta- 
tion function is suppressed at large z above the position of 
the peak. This is in sharp contrast to the case of modified 
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Fig. 4. The nuclear size, Ra, dependence of charm quark 
energy loss for different values of and xb- 
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Fig. 5. Modification factor for the charm quark fragmentation 
function in a nucleus. The value xa ~ 0.05 corresponds to a 
nucleus with a radius Ra = 4.25 fm. 



light quark fragmentation functions which deviate from 
the vacuum form in a very large range of z. 



3 Non-Abelian feature of jet quenching 

Another non-Abelian feature of the parton energy loss is 
its dependence on the color representation of the propa- 
gating parton. The energy loss for a gluon is 9/4 times 
larger than a quark. One can investigate the consequences 
of this non-Abelian feature in the fiavor dependence of the 
high-pr hadron suppression In the mean time, we can 
also study the effect of the non-Abelian parton energy loss 
on the energy dependence of the inclusive hadron spectra 
suppression (23 . One can exploit the well-known feature 
of the initial parton distributions in nucleons (or nuclei) 
that quarks dominate at large fractional momentum (x) 
while gluons dominate at small x. Jet or large pt hadron 
production as a result of hard scatterings of initial par- 
tons will be dominated by quarks at large xt = "^Pt/Vs 
and by gluons at small xt ■ Since gluons lose 9 /4 more en- 
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energy loss g/q=9/4 ■ 
energy loss g/q=1 ■ 
short lifetime t=0-2 fm/c ■ 



Fig. 6. Nuclear modification factor RauAu for neutral pions 
as function of collision energy at fixed pt = 6 GeV in most 
central collisions (with centrality 10%). Here we compare the 
QCD energy loss and a non-QCD one where the energy loss is 
identical for quarks and giuons. 



ergy than quarks, the energy dependence of the large (and 
fixed) pt hadron spectra suppression due to parton energy 
loss should reflect the transition from quark-dominated jet 
production at low energy to gluon-dominated jet produc- 
tion at high energy. Such a unique energy dependence of 
the high-pT hadron suppression can be tested by combin- 
ing i/s = 200 AGeV data with lower energy data or future 
data from LHC experiments. 

To study the sensitivity of hadron spectra suppression 
to the non-Abelian parton energy loss, we calculate the 
single hadron spectra and the suppression factor with two 
different parton energy loss: one for the QCD case where 
the energy loss for a gluon is 9/4 times as large as that 
for a quark, i.e. AEg/AEq — 9/4; the other is for a non- 
QCD case where the energy loss is chosen to be the same 
for both giuons and quarks. Similarly, the average number 
of inelastic scatterings obeys (^)g/(^)g = 9/4 in the 
QCD case. For the non-QCD case we are considering, the 
above ratio is set to one. In order to demonstrate the differ- 
ence between QCD and non-QCD energy loss, we compute 
the Raa for neutral pions at fixed pt = 6 GeV in central 
Au + Au collisions as a function of y/s from 20 AGeV to 
5500 AGeV. Shown in Fig.|Slare the calculated results with 
both the QCD and non-QCD energy loss. In the calcula- 
tion, we fix the initial gluon number density and quarks' 
mean-free path to fit the overall hadron suppression in the 
most central Au -f Au collisions at ^ — 200 GeV. For any 
other energy and centralities, we simply assume the initial 
gluon number is proportional to the final measured total 
charged hadron multiplicity per unit rapidity. One can see 
that due to the dominant gluon bremsstrahlung or gluon 
energy loss at high energy the Raa for the QCD case is 
more suppressed than the non-QCD case where the gluon 
energy loss is assumed to take an equal role as the quark. 

Another interesting feature with the energy depen- 
dence of Raa is the change of slope around = 1300 
GeV. The rapid decrease of Raa at = 20-1300 GeV is 



mainly due to increased initial gluon density and also the 
change of pt slope of jet production cross section with ^/s. 
As the energy loss increases, more jets produced inside the 
overlapped region are completely suppressed. Only those 
that are produced within an out-layer in the overlapped 
region will survive. This will be like surface emission with 
a finite depth. The suppression factor Raa will then be 
determined by the width of the out-layer which is just the 
averaged mean-free-path (A). As a consequence, Raa will 
then have much weaker ^/s dependence. This effect was 
also seen by calculations in Ref. | 28| . 



4 Modified dihadron fragmentation function 

In addition to single inclusive hadron spectra from jet frag- 
mentation, multiple hadron correlations have proven to be 
a useful measurement for characterizing the modification 
of the jet structure in hot medium. For example, while 
the back-side two-hadron correlation is completely sup- 
pressed, the same-side correlation is observed to change 
little when the transverse momentum of the secondary 
hadron is large |21j . Such same-side correlation is essen- 
tially given by dihadron fragmentation functions. It is im- 
portant to investigate why such a dihadron fragmentation 
function changes little in the kinematic region of the ex- 
periments and whether it will be significantly modified 
when the secondary hadron is soft. 

One can define the dihadron fragmentation functions 
in terms of the overlapping matrix between parton field 
operators and the final hadron states, similarly as the 
single hadron fragmentation functions. In the light-cone 
gauge, the dihadron fragmentation function of a quark 
are defined as 



(fphi^ f d'^P 



Tr 



2zh,zhj 2(2^)3 y (27r)4 
X S (^zh - ^) / d'^xe-'P-'= 

s,Ph,,Ph2)(ph2,Phi,s\ij{x)\ 0) 

-^Ph s 



,(3) 



where Zh = Zh^ + Zh^ and pt = Phi + Ph2 ■ Like sin- 
gle hadron fragmentation functions, dihadron fragmenta- 
tion functions also contain non-perturbative physics and 
thus are not calculable in pQCD. However, one can study 
their evolution with the energy scale within pQCD. The 
DGLAP evolution equations for dihadron fragmentation 
functions have been derived recently |29j . For dihadron 
fragmentation function of a quark, it has a form 



dD'^^f^^-{zH,,Zh2,Q^) 

aiogQ2 

2tt 



dz 



dz 



z z 

Zh2 
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2)+ r ^7,.(^K^'^^(^,^,Q^ 



(4) 



that is similar to the DGLAP evolution equations for sin- 
gle hadron jet fragmentation functions. However, there are 
extra contributions in the above equation that are propor- 
tional to convolution of two single hadron fragmentation 
functions. These correspond to independent fragmentation 
of both daughter partons after the parton split in the ra- 
diative correction. Here is the normal parton split- 
ting function and "fqq is the same function without the 
virtual corrections. For numerical solution of the DGLAP 
evolution equations of the dihadron fragmentation, results 
from JETSET [301 at a given scale is used. The evolved 
dihadron fragmentation functions from the DGLAP equa- 
tions at higher scales agree with the Monte Carlo results 
very well |^. It will be useful to compare actual experi- 
mental data when they become available. 

Since the induced bremsstrahlung in medium is sim- 
ilar to that in vacuum, medium modification to the jet 
fragmentation functions should resemble the radiative cor- 
rections in vacuum that lead to the DGLAP evolution 
equations. Therefore, it is not surprising that the medium 
modification to the dihadron fragmentation functions has 
the identical form as the above DGLAP evolution equa- 
tions |31| . These medium modifications depend on the 
same gluon correlation functions as in the modification 
to the single hadron fragmentation functions. Therefore, 
in the numerical calculation of the medium modification 
of dihadron fragmentation functions, there are no addi- 
tional parameters involved. The predicted results are in 
good agreement with HERMES data [311 ■ We find that 
most of the nuclear modification is manifested in the sin- 
gle hadron fragmentation functions. Since dihadron frag- 
mentation functions already contain the information of 
single hadron fragmentation function, the modification to 
the remaining correlated distribution is very small. So 
the normalized correlation D'^^'^'' {zi, z^) / D^^ {zi)D'^-' {Z2) 
has much smaller nuclear modification as compared to 
the single hadron fragmentation functions. This also ex- 
plains why the same-side two-hadron correlation in cen- 
tral heavy-ion collisions remains approximately the same 
as vnp+p collisions while the back-side is completely sup- 
pressed [21]. However, because of trigger bias, one might 
sample different values of zi and Z2 in Au -\- Au and p + p 
collisions. This could lead to apparent change of dihadron 
correlation |5T| . 



5 Jet fragmentation and quark recombination 

During the propagation and interaction inside a decon- 
fined hot partonic medium, a fast parton has not only 
induced gluon radiation but also induced absorption of 
the surrounding thermal gluons. This leads to a stronger 
energy dependence of the net energy loss for an interme- 
diate energy parton |32| . In principle, one can consider 
such processes of detailed balance as parton recombina- 
tion and they can continue until the hadronization of the 



bulk partonic matter. Eventually, during the hadroniza- 
tion, partons from the jet can combine with those from 
the medium to form final hadrons. Indeed, there exists al- 
ready some evidence for the parton recombination in the 
experimental data on the final hadron spectra in heavy- 
ion collisions at RHIC. At intermediate = 2 — 4 GeV/c, 
the suppression of baryons due to jet quenching is signifi- 
cantly smaller than mesons leading to a baryon to meson 
ratio larger than 1, about a factor of 5 increase over the 
value in p + p collisions |33| . On the other hand, the az- 
imuthal anisotropy of the baryon spectra is higher than 
that of meson spectra. Such a flavor dependence of the 
nuclear modification of the hadron spectra and their az- 
imuthal anisotropy is not consistent with a picture of pure 
parton energy loss. The most striking revelation of the un- 
derlying mechanism comes from the empirical observation 
of the scaling behavior between the azimuthal anisotropy 
of baryons and mesons [SJ, {pt/2)/2 = ?;^(pt/3)/3, 
inspired by a schematic model of constituent quark recom- 
bination of hadron production. 

Many quark recombination models "35"36' are success- 
ful in describing the observed flavor dependence of the nu- 
clear modification of hadron spectra at intermediate px- 
These models generally have three different contributions 
to the final hadron spectra. They include recombination 
of the thermal quarks in the bulk matter into hadrons 
which dominate low pt spectra and recombination be- 
tween thermal quarks and quarks from high transverse 
momentum jets that are responsible for intermediate pT 
hadrons. They all assume a thermal distribution for the 
medium quarks and employ constituent quark model for 
the hadron wavefunctions which determine the recombi- 
nation probabilities. However, current models differ in the 
determination of the constituent quark distributions from 
high Pt jets and there exist ambiguities in the connec- 
tion between partons from pQCD hard processes and con- 
stituent quarks that form the final hadrons. Perhaps the 
most consistent treatment of the problem is the model 
by Hwa and Yang |35j . In this model, quark recombina- 
tion processes are traced back to parton fragmentation 
processes in vacuum. They consider the initial produced 
hard partons that will evolve into a shower of constituent 
quarks which then recombine to form the final hadrons 
in the parton fragmentation process. The recombination 
of the shower quarks of the parton jets with the medium 
quarks in heavy-ion collisions can be carried out straight- 
forward given both the shower and medium quark dis- 
tributions. Since this is a phenomenological model, the 
nuclear modification of the jet shower quark distributions 
and their QCD evolution cannot be calculated in the Hwa- 
Yang model. The model has to rely on fitting to the ex- 
perimentally measured hadron spectra to obtain the cor- 
responding nuclear modified shower quark distributions 
for each centrality of heavy-ion collisions and correlation 
between shower quarks are completely neglected. 

We have made a first attempt to derive the quark re- 
combination model for jet fragmentation functions from 
the field theoretical formulation and the constituent quark 
model of hadrons [SZ| ■ Within the constituent quark model, 
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we consider the parton fragmentation as a two stage pro- 
cess. The initial parton first evolves into a shower of con- 
stituent quarks that subsequently will combine with each 
other to form the final hadrons. Since constituent quarks 
are non-perturbative objects in QCD similarly as hadrons, 
the conversion of hard partons into showers of constituent 
quarks is not calculable in pQCD. However, we can de- 
fine constituent quark distributions in a jet as overlapping 
matrices of the parton field operator and the constituent 
quark states, similarly as the definition of the hadron frag- 
mentation functions. 

Given a hadron's (a meson for example) wavefunction 
in the constituent quark model. 



in spirit to the Hwa-Yang recombination model. Given a 
form of the hadrons' wavefunction in the constituent quark 
model, one can in principle extract constituent quark dis- 
tribution functions from the measured jet fragmentation 
functions. Furthermore, within this framework, one can 
also derive the DGLAP evolution equations for the di- 
quark distribution functions, 



\Ph) 



dxi 



2(27r)3 ^x,{l-x,) 
X |fci^,a::i;-fcij.,l - a;i). 



iph{ki±,xi]~ki^,l - xi) 
(5) 



and neglecting interferences between recombination of quark 
and anti-quark pair with different momentum, one can 
rewrite the single inclusive meson fragmentation function 
as a convolution of the diquark distribution functions and 
the recombination probability. 



K 
p 



X Tr 



+ 



{0\m\S,Ph){ph,SUix)\0) 



Jo ^ Zm 



zi) , 



(6) 



where Ch is a constant representing contributions from 
interference processes. The recombination probability is 
determined by the hadrons' wavefunction. 



Rh{kiA_, — ) 

Zh 



Vh{ki^, — ; -fcij^, 1 - — ) 

Zh Zh 



(7) 



and the double constituent quark distribution function is 
defined as the overlapping matrix between the current 
quark operators and the intermediate constituent quark 
states, 



2ziZ2 J 2(27r)3 7(2 



P Li xf Ph 



■^Tr 



^ ^(0 m)\ S, fci, k2){k2,k,, S \i>ix)\ 0) 

/^Ph S 



(8) 



dQ 



:F^'''Hz,,Z2, 



dz 



27r 7^1+^2 z'^ 



dQ- 



;F^^^'^Z,,Z2,Q') 



z z 

2\ /•! 



dz 



27r 7^1+^2 z'^ 



iMFr^i^, ^, Q2) + j^^(z)F^^'>^ (-i, ^, Q2) (10) 



Note that the above equations are a little different from 
the DGLAP evolution equations for dihadron fragmenta- 
tion functions. There is no contribution from independent 
fragmentation in the diquark distribution. This is because 
the diquark distribution functions defined in the context 
of quark recombination are only for two quarks whose rela- 
tive transverse momentum is limited by the wavefunction. 
Therefore, change of the initial momentum scale does not 
lead to variation of phase space available for the diquark 
from the independent fragmentation in the final states. 
One can show that combining the above evolution equa- 
tions for the double constituent quark distribution func- 
tions with the expression of jet fragmentation function in 
Eq. (|14() , the DGLAP evolution equations for single inclu- 
sive hadron fragmentation functions can be recovered. 

The above reformulation of the jet fragmentation func- 
tions does nothing to simplify the complexity of jet hadroniza- 
tion. However, extending the formalism to finite tempera- 
ture, we can automatically derive the contributions from 
recombination between shower and thermal quarks in ad- 
dition to soft hadron production from recombination of 
thermal quarks and leading hadron from recombination of 
shower quarks. The shower and thermal quark recombina- 
tion involves single quark distribution functions which are 
related to the diquark distributions through sum rules. 
Therefore, one can consistently describe three different 
processes within this formalism. Since the single and di- 
quark distribution functions are defined at finite temper- 
ature which are different from the corresponding vacuum 
distributions, one can also consistently take into account 
parton energy loss and detailed balance effect for jet frag- 
mentation inside a thermal medium. 



Here, ph — ki + k2 and Zh — zi+Z2- A is the cutoff for the 
intrinsic transverse momentum of the constituent quarks 
inside a hadron, as provided by the hadron wavefunction. 
The above definition of diquark distribution function in a 
fragmenting parton jet has exactly the same form as the 
dihadron fragmentation functions (29) . One can similarly 
express the parton fragmentation functions for baryons in 
terms of triquark distribution functions. This is similar 



6 Modified jet fragmentation due to quark 
recombination in medium 

One can study the fragmentation of a parton jet in medium 
simply by replacing the vacuum expectation in the S ma- 
trix of the processes or the operator definition of the par- 
ton fragmentation functions by the thermal expectation 
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values, (OIOIO) ^ ((O)), 



m) = 



Tr e-H/T 



(11) 



where, H is the Hamihonian operator of the system and 
T is the temperature. Therefore, the single hadron frag- 
mentation at finite temperature for a quark is defined as 



Tr 



Y,{{m\S,Pk){Pk,S\yj{x))) 



,(12) 



where ph and p are the four-momentum of the hadron 
and the initial parton, respectively. After making all possi- 
ble contraction between the final constituent quark states 
with the thermal intermediate states, one can obtain three 
distinctive contributions to the above fragmentation func- 
tion in medium. 



D^^iz,,p+) = D'^(^^\z,,p+)+D^(''^\z,,p+) 
+ D^^(^^Hz,,p+). 

The first contribution, 



(13) 



D 



h(SS) 



dzi 



RhiO 



X F;i'^^{zi,Zh 



^) 

Zh 

zi,P'^) , 



(14) 



normally referred to as " shower-shower" quark recom- 
bination comes from recombination of constituent quarks 
from within the parton jet. It has exactly the same form 
as the parton fragmentation functions in vacuum [Ea. H14|) ] 
in the framework of quark recombination, except that the 
diquark distribution function F^^^^^ {zi, Z2) are now also 
modified by the medium. Its definition is similar to that in 
vacuum in Eq. JHJ but the vacuum expectation is replaced 
by thermal average. These modified diquark distribution 
functions should in principle contain effects of multiple 
scattering, induced gluon radiation and absorption, in the 
same way as the modification of hadron fragmentation 
functions in a thermal medium j38j . 

The second term in the modified fragmentation func- 
tion, 



fq{q±,ZqP+)F^{zh 



dz„ 



(fq±_ Rh{q±,Zq/zh) 



Zh J 2(27r)3 (1 - Zq/zhf 

Zq) + fq{q±,ZqP+)F^{zh - Zg) (15) 



is from recombination between a constituent quark (anti- 
quark) from the parton jet and an anti-quark (quark) from 
the medium. This is often referred to as "thermal-shower" 
quark recombination. Here, fq and fq are thermal quark 
distributions, and Flj{z) and F^{z) are single constituent 
quark or anti-quark distributions of the fragmenting par- 
ton jet in a thermal medium defined similarly as in the 



vacuum, except that the vacuum expectation values are 
replaced again by thermal averaged expectation. They 
should be different from the corresponding quark distri- 
butions in vacuum because of multiple scattering, induced 
gluon bremsstrahlung and parton absorption. 
The final term in Eq. ((T^ . 



D'^q^^^\zh,P+) = Vp^ 



(27r)S 



dZq 



(27r)3 



fqiq^, ZqP'^)fq{Phl^ - 9i , {z,, - Zq)p+)Rh{q^, (16) 

Zh 

comes from recombination of two thermal constituent quarks. 
Here V is the total volume of the whole thermal system. 
Since hadron production from the thermal quark recom- 
bination is not correlated with the parton jet and its frag- 
mentation, the above expression is a little bit unnatural. 
One should be able to rewrite it in terms of the invariant 
hadron spectrum from thermal quark recombination. 



(2^) 



dp'^d^Ph 



= V dxn 



d^ 



q± 



„ .'h± Jo J 2(27r) 

^fq{Ph± - q±, (1 - Xq)p+)Rh{q^,Xq) . 



fq{q- 



'±:^qPh) 



(17) 



which is not correlated and therefore do not depend on 
the parton jet fragmentation. The above expression also 
coincides with results from other recombination models 

m 



7 Conclusions 

In this talk, I reviewed some new developments in the 
study of modified single and dihadron fragmentation func- 
tions in dense medium and their applications to heavy-ion 
collisions. All of them will help to provide further test of 
the picture of parton energy loss and jet quenching and en- 
able more detailed characterization of the dense medium 
those jets probe. The mass dependence of the gluon forma- 
tion time from the heavy quark leads to a unique change 
of the medium size dependence of the heavy quark energy 
loss, from linear to quadratic, when the initial quark en- 
ergy and the momentum scale are varied. The so-called 
"dead-cone" effect, also caused by the heavy quark mass, 
in addition reduces the total heavy quark energy loss. The 
form of heavy quark fragmentation function into heavy 
quark mesons in vacuum, which is peaked at large frac- 
tional momentum z, leads to a medium modification that 
is different for light hadrons from massless partons. Since 
one can identify heavy quark mesons, one can use them 
to tag heavy quark propagation and study the difference 
between quark and gluon energy loss |39| . One can also 
use the energy dependence of the suppression of single 
inclusive hadron spectra to test the difference in quark 
and gluon energy loss due to the non-Abelian gauge in- 
teraction. For fixed px the fraction of initial jets changes 
from quark to gluon dominated partons and the different 
energy loss of quarks and gluons will result in a unique 
energy dependence. 



8 



Xin-Nian Wang: QGP and Modified Jet Fragmentation 



Finally, I also discussed the formulation of fragmen- 
tation functions in a quark recombination picture within 
a constituent quark model. Given the hadron's wavefunc- 
tion in the constituent quark model and neglecting inter- 
ference effects, we have shown that hadron fragmentation 
functions can be expressed as the convolution of the re- 
combination probability (given by the hadron's wavefunc- 
tion) and the constituent (or shower) quark distribution 
of the jet. The constituent quark distributions are defined 
as the overlapping matrices between parton fields and the 
final constituent quark states, just like hadron fragmen- 
tation functions as overlapping matrices between parton 
fields and final hadrons. We have derived the DGLAP evo- 
lution equations for the quark distribution functions. We 
then extended the formalism to include the medium effect 
within the framework of field theory at finite temperature. 
One naturally arrived at three distinctive contributions 
from recombination between shower constituent quarks, 
shower-thermal and thermal-thermal quarks, as have been 
proposed by previous recombination models. 

Since parton energy loss is only sensitive to the ini- 
tial color charge (or gluon) density of the medium, the 
observed jet quenching points us to an enormously high 
initial gluon density created in the central Au + Au colli- 
sions at RHIC. However, parton energy loss does not dis- 
tinguish confined and deconfined matter. Quark recom- 
bination between jet shower quarks and thermal quarks 
on the other hand is a process of quark interaction that 
crosses hadronic boundary. If proven, the combined sig- 
nal of quark recombination and thermalization will lead 
to unambiguous conclusion of deconfinement of the pro- 
duced dense matter. 
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tor, Office of Energy Research, Office of High Energy and 
Nuclear Physics, Division of Nuclear Physics, and by the 
Office of Basic Energy Science, Division of Nuclear Sci- 
ence, of the U.S. Department of Energy under Contract 
No. DE-AC03-76SF00098. 
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